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established.  It  consists  of  combining  the  observed  phases  linearly  with  weights 
depending  upon  the  element  locations.  It  is  shown  that  the  presence  of  thermal 
noise,  for  sufficiently  high  signal  to  noise  ratio,  does  not  change  the  rtructure 
of  the  estimator.  Comparison  with  conventional  multiple  interferometric 
techniques  indicates  the  superiority  of  the  proposed  scheme.  Finally,  a 
limited  numerical  study  on  a small  linear  array  vertically  located  on  a reflect- 
ing terrain  is  performed.  Although  in  such  a situation  the  scheme  proposed  is 
not  the  theoretical  optimum,  it  leads  to  errors  that,  for  most  directions  of 
the  target,  are  smaller  than  those  found  for  the  same  array  when  using  con- 
ventional multiple  interferometer  techniques. 
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Maximum  Theoretical  Angular  Accuracy  Of 
Planar  and  Linear  Arrays  Of  Sensors 


1.  INTRODUCTION 

The  objectives  of  this  paper  are: 

To  establish  an  upper  bound  on  the  theoretical  accuracy  in  the  angular 
location  of  a radiating  object  obtainable  with  an  array  of  identical  sensors  subject 
to  instrumental  errors,  in  the  absence  of  angular  interference; 

The  determination  of  an  optimum  processing  scheme,  which,  under  cer- 
tain assumptions,  meets  the  bound; 

A limited  numerical  investigation  of  a particular  array  in  the  reasonably 

realistic  situation  for  which  terrain  specular  reflection  is  present. 

The  question  of  the  errors  in  determining  the  angular  coordinates  of  an  object 

located  in  the  far  rone  of  a planar  or  linear  array  of  antennas  has  been  studied, 

1-3 

among  other  applications,  for  radar  antennas  and  for  multiple  interferometers 
for  landing  systems.4  In  the  first  area  of  applications  mentioned,  the  array 
(Received  for  publication  29  April  1976) 

1.  Rondinelli,  L.  A.  (1959)  Effects  of  random  errors  on  the  performance  of 

antenna  arrays  of  many  elements,  1959  IRE  Nat.  Conv.  Record,  pt  1, 
pp  174-189.  

2.  Allen,  J.  L.  (1963)  The  Theory  of  Array  Antennas.  MiT  Lincoln  Laboratory 

Tech  Report  No.  323,  pp  69-7d. 

3.  Carver.  Keith  R.,  Cooper,  W.K.,  and  Stutsman,  W.L.  ( 1973) Beam  pointing 

errors  of  planar  phased  arrays,  IEEE  Transaction  on  Antennas  and  Propaga- 
tion,  A P-21.  No.  2:199-202. 

4.  Danville,  A.  R.  and  Moore,  S.  R.  (1975)  Design  method  for  interferometer 

arrays.  Proceedings  of  the  IEE  (British).  122,  No.  6:605-#08.  * 
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element  outputs  are  linearly  combined  in  the  feed  structure  in  order  to  provide  a 
small  set  of  "observable"  voltages  (the  outputs  of  the  sum  and  difference  ports  of 
a monopulse  antenna  system)  from  which  the  unknown  target  angular  coordinates 
are  to  be  extracted.  In  the  second  area  of  application  the  angular  position  is  ob- 
tained by  averaging  phase  measurements  pertinent  to  interferometric  baselines 
of  different  lengths. 

The  two  techniques  are,  of  course,  totally  different  and  are  applied  also  in 
totally  different  circumstances.  The  first  technique  uses  a "conventional"  and 
hence  linear  antenna.  The  second  instead  embodies  processing  schemes  which 
are  highly  nonlinear.  Also  in  most  cases  the  second  technique  is  applied  to  lo- 
cate only  cooperative  targets,  (as  in  microwave  landing  systems!  having  a trans- 
ponder on  board.  Thus  the  density  of  the  received  power  incident  on  the  array  can 
be  orders  of  magnitude  greater  than  for  the  radar  case,  thereby  making  considera- 
tions of  antenna  gain  totally  different  in  the  two  cases.  In  this  work  the  question 
of  the  angular  accuracy  of  a system  of  sensors  having  equal  radiation  patterns  is 
reexamined  from  a different  and  more  general  viewpoint.  Rather  than  analysing 
a particular  processing  scheme  assumed  a priori,  an  "unstructured"  approach  is 
adopted  based  on  statistical  inference  considerations.  An  estimation  criterion 
rather  than  a processing  scheme  is  established  a priori  for  the  extraction  of  the 
angular  parameters  from  the  phases  of  the  complex  voltages  observed  at  the  ele- 
men*  ports  (the  amplitudes  carrying  r.o  angular  information).  The  processing 
scheme  is  instead  obtained  from  the  set  of  equations  that  are  the  analytical  ex- 
pression of  the  criterion.  The  maximum  likelihood  criterion,  with  its  well-known 
desirable  statistical  properties,  is  the  one  chosen.  If  the  probability  densities  of 
the  phase  errors  of  the  element  voltages  are  normal,  the  estimator,  that  is,  the 
processing  scheme  of  the  observables,  consists  of  linearly  combining  the  phases 
of  the  element  outputs  with  weights  depending  upon  the  element  coordinates,  re- 
ferred to  the  geometrical  "center  of  mass"  of  the  array  (Section  3).  For  normal 
density  of  phase  errors  the  estimator  is  "optimal",  in  the  sense  of  having  the 
smallest  possible  ms  estimation  errors  among  all  possible  estimators  that  are 
"unbiased",  that  is,  without  systematic  errors.  It  is,  in  fact,  shown  that  under 
the  hypothesis  of  gaussian  probability  densities  for  uncorrelated  errors  in  the 
measurements  of  element  phases,  the  rms  estimation  errors  meet  a Cramer- Rao 
bound  derived  in  Appendix  C.  Thus  the  proposed  algorithm  provides  the  maxi- 
mum possible  angular  accuracy  (for  unbiased  estimation),  which  cannot  be  im- 
proved no  matter  how  sophisticated  the  design  of  the  processing  scheme. 


5.  Van  Trees,  H.  L.  (196b)  Detection  Estimation  and  Linear  Modulation  Theory. 

John  Wiley  and  Sons,  Chap.  1. 

6.  Van  Trees,  H.  I.,  (1968)  op.  cit. , Chap.  2. 
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In  this  paper  planar  and  linear  arrays,  with  elements  in  general  non-equally 
spaced,  are  considered.  As  a numerical  application,  the  performance  of  a linear 
array,  having  the  elements  spaced  in  geometric  progression  (with  the  minimum 
distance  between  elements  small  enough  to  eliminate  ambiguity  in  the  phase  read- 
ings) with  a total  length  of  24  wavelengths,  is  discussed.  The  accuracy  found  is 
substantially  greater  than  that  obtained  for  the  same  geometry  and  phase  errors 
by  using  multiple  interferometer  techniques  (Section  4).  For  the  same  array  ver- 
tically located  on  a reflecting  terrain,  the  degradation  of  performance  due  to  multi- 
path  is  also  examined  (Section  5).  For  this  situation  the  estimation  procedure  pro- 
posed is,  of  course,  not  the  theoretical  optimum.  In  fact,  not  only  phases  but  also 
relative  amplitudes  of  element  output  voltages  should  be  used  because  they  also 
carry  angular  information.  Also  in  this  ease  the  phase  "errors"  due  to  multipath 
are  highly  correlated  and  difficult  to  characterize  statistically.  Nevertheless,  for 
the  particular  case  considered,  the  performance  of  the  scheme  discussed  here  is 
substantially  better  than  that  obtainable  through  a conventional  multiple  interferom- 
eter technique.  This  leads  to  the  conjecture,  to  be  verified  through  analysis  and/or 
simulation,  that  the  scheme  is  in  general  superior  to  the  multiple  interferometer 
technique  also  in  the  realistic  ease  of  presence  of  reflection. 

Finally,  it  is  shown  that  the  inclusion  in  the  model  of  thermal  noise  (assumed 
absent  when  obtaining  the  above  results)  does  not  affect  the  structure  of  the  esti- 
mator (but,  of  course,  affects  the  accuracy)  (Section  6). 


2.  STATEMENT  OF  THE  PROPLEM  AND  NOTATION 

The  discussion  is  conducted  in'terms  of  planar  arrays.  The  results  for  the 
linear  simpler  case  are  then  obtained  from  those  pertinent  to  the  planar  case 
through  appropriate  simplifications.  Our  main  objective  is  to  assess  the  effects  of 
instrumentation  errors,  that  is  the  errors  affecting  the  voltages  at  the  element  ports. 
Thus  (at  least  in  the  first  part  of  this  study)  thermal  noise,  introduced  in  the  voltage 
measurements,  that  is  noise  at  the  element  receivers,  will  be  considered  negligible 
with  respect  to  the  incoming  signal.  The  elements  of  the  array,  or  rather  their 
radiation  patterns,  are  assumed  identical  and  identically  oriented.  Since  (at  least 
for  the  time  being)  we  restrict  our  considerations  to  the  angular  location  of  a single 
target  in  the  absence  of  angular  interferences,  no  angular  information  is  contained 
in  the  amplitudes  of  the  element  output  voltages.  Thus  only  the  set  of  the  phaues  of 
the  element  voltages  constitutes  the  observables  to  be  used  for  the  statistical  ex- 
traction of  the  angular  coordinates. 

I.et  Xj,  >'j  be  the  coordinates  of  the  ith  sensor  in  the  x,  y plane  of  he  array, 
the  total  number  of  elements  being  N.  It  proves  convenient  to  introduce  the  coor- 
dinates of  the  geometrical  "Center  of  Mass"  of  the  antenna  elements: 
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x = a 

£ x.  . 
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(1) 
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The  following  geometrical  quantities  are  also  introduced  for  formal  simplification 
of  our  results: 
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(x.  - x)  , 
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Let  a,  0 be  the  unknown  cosines  of  the  direction  of  incidence  with  respect  to  the 
x,  y axis.  If  k = 2«/X  i8  the  free-space  propagation  constant,  the  propagation  vec- 
tor of  the  incident  wave  is 

k = k|a  x + 0 y + >/ 1 - a ' - 0 ^ i 
where  x.  y.  i are  unit  ve  tors  in  the  directions  of  the  rectangular  axes. 

The  phase  ^ observed  at  the  output  of  the  ith  element  is  expressed  at  the  sum 
of  three  terms: 


vi  * k(“  xi  + 0 y^  * y • 


(6) 


The  first  term  depends  upon  the  element  position  and  contains  the  unknown  param- 
eters. The  second  term  p is  a reference  phase,  interpreted  t s the  phase  that  one 
would  measure  at  the  output  of  an  element  assumed  located  at  the  origin  of  th>»  coor- 
dinate system  in  the  absence  of  measurement  errors.  The  quantities  6.  denote  the 
instrumental  errors  in  the  phase  measurement.  They  may  depend  upon  mechanical 
tolerances  of  the  elements,  quantization  in  the  phase  measuring  device,  and  so 
forth,  and  are  constant  in  time.  The  important  point  is  that  the  flj's  are  unknown, 
except  for  certain  statistical  properties,  and  can  be  treated  as  random  variables. 

We  will  make  several  reasonable  assumptions  about  the  ftj's.  The  probability 
density  is  the  same  for  all  elements: 
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Pi<6i>  = P<6j)  . 

that  is,  independent  of  i.  The  6j's  have  zero  mean: 
Effij]  = / pCfij)  d 6t  = 0 . 


where  £(  ] represents  the  statistical  average  operator.  Also  errors  in  different 

elements  are  statistically  independent,  that  is,  with  self-explanatory  notation: 

P^fl j,  •••  “ pi6j)  pifig^  •••  P^j^^  i (8) 


so  that: 

£[«,«„]  =»,2«lk.  (9) 

where  is  the  rms  value  of  the  phase  error  and  is  Kronecker's  delta,  equal  to 
1 for  i = k,  equal  to  zero  for  i 4 k.  Because  of  (6),  Eq.  (8)  can  be  written  ns  a con- 
ditional probability  density  for  the  set  of  phases  j j for  given  values  of  the  param- 
eters or,  0 , fi: 

Pt  (fy.  *2  •••  % I “•  0.  M)  * jl  p£*>.  - n - k(a  x.  + 0 y.)j  . CO) 

If  the  function  p(- ) in  Eq.  (8)  is  known,  Eq.  (10)  summarizes  all  the  statistical  in- 
formation available  on  or,  0 , given  the  set  of  observables  { V- 1 • In  the  particular 
case  of  normally  distributed  phase  errors. 


p(6  )=_J exp  ( — 

\ 2nb  J 


We  want  to  establish  a statistically  satisfactory  estimation  algorithm  for  the  ex- 
traction from  the  set  of  observables  of  the  unknown  parameters  or,  0. 


r’— 


3.  MAXIMUM  LIKELIHOOD  ME  rHOD  OF  PHASE  FRONT 
SLOPE  EVALUATION 

3.1  Theory 

Under  the  idealized  assumption  of  absence  of  directional  or  diffused  interfer- 
ence the  phase  front  of  the  wave  incident  on  the  array  is  planar.  Since  all  angular 
information  is  contained  in  the  slope  of  the  phase  front,  an  intuitively  appealing 
procedure  for  the  extraction  of  a,  0 , consists  of  measuring  the  phases  of  the  ele- 
ment voltages  then  combining  them  linearly  to  extract  the  components  of  the  phase 
slope.  If  the  weights  used  in  the  linear  combination  are  properly  chosen,  this  pro- 
cedure implements  the  Maximum  Likelihood  Estimation  (MLE)  criterion  under  the 
assumption  (11)  of  gaussian  probability  density  for  phase  errors.  Furthermore, 
under  the  same  assumption,  it  yields  an  unbiased  estimator  having  the  minimum 
possible  rms  errors. 

We  recall  that  the  MLE  consists  in  assuming  that,  for  a given  set  of  observ- 
ables, the  most  likely  values  of  the  unknown  parameters  are  those  for  which  the 
probability  of  errors  is  maximum.  Hence,  following  the  standard  procedure,  we 
take  the  derivatives  of  the  logarithm  of  the  conditional  probability  densities  (10), 

6 * AAA 

and  equate  them  to  zero.  The  estimates  of  a,  0,  n,  are  those  values  a,  0,  fi 
for  which  the  equation 

= 0 

(12) 

a=a.  0 = 0. 

holds,  together  with  the  two  parallel  equations,  obtained  by  replacing  3a  in  (12)  with 
30  and3(i,  respectively.  Assume  normal  densities  for  the  fl^'s.  Thus  from  (10) 
and  (12)  we  obtain: 
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1*1 
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- kdtXj 

- K0  y,  - 

£)  x.  = 0 , 

(13) 
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(*i 
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* k0  yt  ’ 

y4  » 0 - 

(14) 

1*1 
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i*i 

% 

- kax. 

i 

* k0  y.  - 

= 0 . 

(15) 

By  eliminating  in  (13-15)  the  estimate  /2  of  the  parameter  fi  (of  no  interest)  the 
explicit  forms  of  the  estimators  a and  0 are  obtained  (Appendix  A): 
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In  Appendix  B it  is  shown  that 

E[a]  = a,  E[0]  = 0, 


(18) 


that  is,  the  estimates  06-17)  are  unbiased.  In  Appendix  B the  quadratic  errors  of 
estimations  are  shown  to  be: 


2 "6 
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k M M - M 
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(20) 
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In  Appendix  C it  is  shown  that  under  the  hypothesis  of  normal  densities  tor  the 
phase  measurement  errors  the  rms  estimation  errors  (19-20)  are  the  minimum 
possible  for  unbiased  estimators.  This  means  that  the  algorithm  (16-17)  is  optimal 
for  yielding  the  maximum  possible  angular  accuracy  in  the  absence  of  systematic 
errors. 

* A 

The  cross  correlation  of  the  estimators  a and  0 is  (Appendix  B): 
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From  (19-21)  the  joint  probability  density  for  flf  and  0 is  easily  established: 


p(ot.  0) 


2*  a , 


J M M - M2  exp  [ - [m  (a -ci 
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(a - a)<0  -0)]|  . 


(22) 


By  invoking  the  central  limit  theorem  Eq.  (22)  is  approximately  valid  even  if  the 
p (6j)  are  non-gaussian. 
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3.2  Remark* 

The  phase  reference  is  immaterial,  provided  of  course  it  is  the  same  for  all 
elements.  In  fact,  by  using  the  obvious  identities  (A4)  of  Appendix  A it  is  found 

A A 

that  a and  0 remain  unchanged  ifin(16-17)  every  is  replaced  by: 


m.  - f 
1 1 c 


with  a constant  phase  (independent  of  the  index  i).  For  example  the  phases  can 
be  measured  with  respect  to  the  output  of  an  auxiliary  sensor  whose  position  with 
respect  to  the  array  elements  is  immaterial  and  need  not  even  be  known.  Alterna- 
tively and  more  conveniently  the  phases  can  be  measured  with  respect  to  a particu- 
lar array  element,  say  the  one  indexed  as  1.  This  leads  to  an  estimator  whose 
expression,  while  formally  different,  is  completely  equivalent  to  (16-17)  because 
its  expression  is  obtained  from  (16-17)  through  analytic  manipulation.  In  fact  the 
following  identity  holds  as  a consequence  of  (A4): 

N 

2 (C,  [<x;  - x)  My  - (yt  - y)  Mx  j 

N N 

= *i  [(xi  ‘ x)  My  ■ (yi  ■ y>  Mx]  -£2  *l[  Ui  - *>  My  • (yi  • y)  Mx]  ‘ 


and  thus  (16)  can  be  written: 


£ <^i-^)[(Xi-x)My-(yi-y)Mx] 


MM-  M 
x y xy 


and  a parallel  expression  can  be  established  for  . Since  (23)  is  identical  to  (16) 
all  the  statistical  properties  established  for  (16)  evidently  hold  for  (23)  also.  Of 
course  in  this  case,  unlike  the  scheme  using  an  additional  element  for  phase  refer- 
ence, the  position  of  the  element  index  by  1 plays  a fundamental  role  in  the  algor- 
ithm (23)  because  it  appears  in  the  quantities  (1-5). 

Since  the  phases  arc  measured  modulo  2ir  there  i1'  an  inherent  ambiguity  in  the 
phase  measurement.  This  can  be  eliminated,  however,  by  choosing  the  spacings  of 
some  elements  small  enough  to  yield  an  unambiguous  phase  difference  between  con- 

4 

tiguous  elements.  The  technique  is  in  principle  well  known  and  a detailed  study  of 
this  point  is  outside  the  scope  of  this  paper.  However  a brief  discussion  will  be 
made  in  Appendix  D in  connection  with  the  numerical  example  of  Section  4. 


i 


It  is  reasonable  to  assume  that  in  practical  cases  the  array  will  have  at  least 
one  axis  of  symmetry,  x or  y.  This  assumption  makes 


M 


xy 


o 


A * 

[see  (5)].  The  two  estimates  a and  0 are  statistically  independent  fas  shown  by 
(221  and  thus  (16)  takes  le  form: 

N 


L ip  (x.  - 5) 
1 1=1  1 1 
k — *r — 


(24) 


A similar  form  holds  for  0 (with  y replacing  x).  The  rms  estimation  error  is: 


2 

2 1 

= ~T  W~  ' 
k x 


(25) 


2 

with  a parallel  expression  for  a ^ . Expressions  (24-25)  are  of  course  also  those 
holding  for  linear  arrays.  The  one-dimensional  version  of  the  modified  algorithm 
(23)  is 


N 

T/  iip.  - »>.)  (x.  - x) 

» t i-  2 1 

= E • (26) 


It  is  clear  from  (19-20)  or  (25-26)  that  to  increase  accuracy  for  a given  number 
of  elements  it  is  convenient  to  have  the  elements  as  sparse  as  possible,  compatible 
with  the  need  of  avoiding  phase  ambiguities,  in  order  to  increase  the  "moments  of 
inertia"  M . M , directly  affecting  the  accuracy  of  the  estimates  of  the  angular 

X y 

location. 


4.  LINEAR  ARRAYS.  COMPARISON  WITH  ACCURACY 
OF  MULTIPLE  INTERFEROMETERS 

The  processing  scheme  discussed  in  the  preceding  sections  provides  consider- 
ably more  accuracy,  that  is,  smaller  rms  angular  estimation  error,  than  does  the 
"conventional"  multiple  interferometer  technique.  The  comparison  is  made  by  con- 
sidering the  same  linear  array  (that  is  with  the  same  element  locations  ar.d  measure- 
ment errors)  using  two  different  processing  schemes:  the  one  proposeJ  here  and 
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that  used  in  the  "conventional"  interferometric  technique.  In  the  latter  the  angular 
information  is  extracted  through  the  algorithm  described  in  Ref.  4,  which  with  our 
notation  takes  the  form: 


N 


“m  ' E TT 


E ( ii> . - *,) 

1 i=2  1 1 


(27) 


E (x.  - x.) 
1=2  1 1 


The  estimator  (27)  is  unbiased.  Its  variance  is  found  to  be: 
2 


2 

°m  = 


N(N  - 1) 

TT 

^ 


2 * 


(28) 


Expression  (28)  is  always  greater  than  (25).  The  result  can  be  established  simply 
by  recalling  that,  for  gaussian  probability  densities  of  the  phase  measurements, 
the  rms  error  in  angular  location  associated  with  the  estimator  (24)  or  (26)  is  the 
minimum  possible  for  the  given  array  geometry,  because  it  meets  the  bound  on  rms 
estimation  error  in  Appendix  C.  However,  it  can  be  found  more  directly  by  using 
simple  arithmetic  inequality. 

In  order  to  gain  a better  appreciation  of  the  superiority  of  the  estimator  pro- 
posed in  comparison  with  a multiple  baseline  interferometer  scheme,  consider  the 
following  particular  case. 

Example 

Suppose  that  it  is  known  a priori  that  the  airection  of  the  object  is  within  a 
sector  of  ± 15*  from  the  plane  perpendicular  to  the  array  axis.  In  order  to  avoid 
ambiguities  we  will  chose  the  minimum  spacing  among  the  elements  equal  to  0.  75  ) . 

The  distances  between  first  and  successive  elements  increase  in  geometric 
progression.  In  this  way  all  the  phase  ambiguities  are  eliminated  through  a simple 
algorithm  discussed  in  Appendix  D.  Thus  if  c is  the  distance  of  the  first  element 
from  the  origin  of  the  x axis,  the  abscissae  of  the  elements  are: 


In  practical  cases  this  may  mean  that  through  an  auxiliary  antenna  system  the  sig- 
nal is  blanked  if  incident  from  outside  the  sector  of  interest,  the  angular  measur- 
ing system  being  aciivited  only  for  targets  within  the  sector. 
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X1  = c 


x2  = c + 0.  75  X , 


Xg  = c + 1.  5 X , 


x^  = c + 3 X , 


x,.  = c + 6 X , 


Xg  = c + 12  X , 


x^  = c + 24  X . 


The  abscissa  of  the  geometric  center  of  mass  is: 


7 

T.  x . 


x = -My — = c + 6. 75  X . 


Its  moment  of  inertia  with  respect  to  the  center  of  mass  is,  from  (3): 


M + 448.  875  X . 
x 


Suppose  that  the  rms  error  in  the  phase  measurements  at  the  element  ports  is: 


Og  = 0.5  radians  - 28.648*  . 


From  (25)  the  rms  error  in  the  object  angular  location  is,  in  or  units: 
2 

1 

TT 


2 "fl 


» (2.)2  N 


1=1 


Put 


a = sin  9 , 


with  9 the  angle  from  the  array  normal.  The  rms  error  of  9 is  maximum  for 
9 = 15*.  the  edge  of  the  sector  of  Interest,  and  is  in  milliradians. 


°9  ~ cos  15 


“ r 103  - 3.89  mrad  . 


(29) 
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Let  us  consider  the  rms  error  <jp  for  the  multiple  interferometer  scheme.  By 
applying  (28)  we  obtain  cq  = 11.3  mrad.  The  rms  error  ?s  substantially  greater 
than  for  the  optimum  scheme. 

If  a different  element  is  chosen  as  a reference  the  rms  error  will  vary,  still 

remaining,  of  course,  greater  than  the  value  (29).  For  example,  assume  as  a 

reference  element  the  one  located  at  x = c +24X.  In  such  a case  the  rms  estima- 
// 

tion  error  becomes  <jp  - 10.6  mrad,  again  greater  than  (29),  as  it  must  be.  We 
recall  that  we  have  considered  the  ideal  situation  in  which  any  angular  interference 
is  absent. 


5.  SPECULAR  REFLECTION  EFFECTS-A  NUMERICAL 
EXPERIMENT 

In  the  presence  of  coherent  interference  like  specular  reflection,  the  proces- 
sing scheme  here  proposed,  although  not  the  theoretical  optimum,  on  the  basis  of 
a limited  numerical  study  seems  to  exhibit  comparatively  good  performance  with 
respect  to  a conventional  multiple  baseline  interferometric  system.  A theoretical 
analysis  seems  impossible  or  at  least  exceedingly  difficult,  mainly  because  phase 
"errors"  induced  by  specular  reflection  are,  of  course,  highly  correlated.  A 
numerical  experiment  made  for  a particular  case  shows  the  superiority  of  the  al- 
gorithm (24)  or  (26)  proposed  here  with  respect  to  the  conventional  multiple  inter- 
ferometer scheme  (27). 

Consider  the  array  geometry  discussed  in  the  previous  section.  Suppose  the 
array  is  vertically  located  on  a reflecting  terrain.  The  distance  of  the  first  and 
lowest  element  from  the  origin  located  on  the  ground  is  c.  If  0 is  the  elevation 
angle  of  the  scattering  object,  and  if  R is  the  magnitude  of  the  reflection  coefficient 
the  output  voltage  at  the  ith  element  is: 

(x.  + c)  sin  Or  j V,  - j^r-  (x.  + c)  sin  6 T 

V.  = e A 1 [l  + Re  1 A 1 J,  (30) 


where  <p j is  the  phase  of  the  reflection  coefficient.  Put 
<p  = Vj  + ^ (c  + x)  sin  8 . 

The  phases  of  the  Vj  are: 

Vj  = -y  (Xj  + c)  sin  0 + tg 


j R sin  [v>  • ’ x>  sin  0 j j 

( 1 ► R cos  f - y-  (Xj  - x)  sin  $1  J 


(31) 


(32) 


Inserting  (32)  into  (24)  or  (26),  one  obtains  for  the  error  in  the  estimation  of  sin  8 


*?  *"•■ ' iww*.?  ■ vmiimmfmmmmm wwjwphw wjuw^wa" 


Ml  UMiMMSI 


N 


£ <xi 


sin  5 - sin  9 = 


1=1 


.(  R sin  [m  - (x.  - x)  Bin  9 1 

‘i>tg"V+ 


R cos  [ifi  - ^ (Xj  - x)  sin  0] 
~?M 

A X 


(33) 


Assume  the  error  is  small.  We  obtain  then: 


N 


_ i ( R sin  [to  - (x.  - x)  sin  0 ] ) 

V (x.  -5)  tg  1 -X  - K- 

| 1 + R cos  (xj-x)  sin  5]  | 


A0  - — 


~2i~ 

X 


(34) 


In  (33),  <fi  is  a phase  angle  depending  upon  the  distance  of  the  array  from  the  re- 
flecting plane,  the  phase  of  the  reflection  coefficient  and  the  element  spacings.  In 
practical  cases  V is  not  known.  Hence  to  assess  the  merit  of  the  processing  scheme, 
it  is  reasonable  to  compute  the  error  as  a function  of  9 for  different  values  of  V and 
a large  R.  In  Figures  1 and  3 the  error  has  been  plotted  for  R = 0.  3 and  values  of 
<(>  equal  to  0,  *,  and  ± w/2.  Figures  2 and  4 show  the  error  in  a multiple  inter- 
ferometer system  with  the  same  array  geometry  and  the  same  values  of  R and 
It  is  apparent  that  for  the  scheme  proposed  here  the  maximum  error  is  less  and 
decreases  more  rapidly  with  increasing  elevation  than  for  conventional  multiple 
interferometer  systems.  Although  this  may  suggest  that  in  general  the  algorithm 
(24)  or  (26)  is  superior  to  (27)  even  in  the  presence  of  specular  reflection,  more 
extensive  computations  are  of  course  necessary  to  corroborate  this  contention. 


Figure  1.  Proposed  Scheme,  Angle  Estimate 
Error,  R * Q.  9,  0 * *,  and  0*0 
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Figure  3.  Proposed  Scheme,  Angle  Estimate 
Error,  R = 0.  9,  6 = y,  0 ’ - j 


6.  THERMAL  NOISE 

We  consider  again  the  situation  of  no  angular  interference.  The  presence  of 
thermal  noise  can  be  accounted  for  by  assuming  that  the  output  voltage  of  the  ith 
sensor  is  a function  of  time  described  by: 

v,(t)  = f A + m.(t)  1 coa(w  t + ♦,)  + n. (t)  sin  lut+  (!>  ) , (35) 

1 L oil  oi 

where  wQ  is  the  center  frequency  and  m.(t)  and  n.(t)  are  assumed  to  be  narrowband 
gaussian  processes.  If  B is  the  bandwidth  (assumed  rectangular),  (35)  can  be 
sampled  at  a sampling  rate  of  s/B.  The  phase  and  quadrature  components  of  the 
samples  are  the  real  and  imaginary  parts  of  the  complex  quantities 


Assume  that,  with  probability  close  to  unity: 
A >>m  (t)  . 

A > >n.(t)  . 
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(38) 


ikJaA-fc*  it’il  % ^ 


I 


Then  (36)  yields 


\ (I) : A“»  SJ  [*1  + i ni  (I)  t"o(l)]j_ 


(39) 


which  will  be  written  concisely  as; 

(JL>  s 


Vl.  * A e*p  [j  (*»ia  )j 


where  the  <p.'s  are  the  phases  of  (39)  and  are  independent  random  variables  whose 

1 C 

probability  densities  are  induced  by  those  of  6.  and  n.(=).  The  former  is  assumed 

1 1 rs 

gaussian  ?nd  given  by  (11).  The  second  is  similarly  given  by: 


P (n.  ) = —1—  exp  (-  n?  / 20  2) 

n is  j2no  \ is  n f 

y n 


(40) 


where  we  have  concisely  written: 
ni  (!)  • 


n. 

is 


Let  M = TB  be  the  total  number  of  samples,  T being  the  observation  time,  The 
probability  of  the  set  of  N x M phases  v>jg»  conditioned  by  the  parameters  a.  0 , fi. 


is: 


Pt  (<PH*  *®12*  * * ^1M’  *21*,,*2M,*N1,,*,9NM 
N TB 


= f=l  J=1  L P6^8  * M - kor  x,  - k 0y,  - p(nl8>  dnia 


or  from  (11)  and  (40),  through  some  manipulation: 


Pt  ^ir ‘”12  •••  o’im*  ^21*  v522,,,p2M' 
N TB 


a.  0,  p) 


fl  II 


1 


1 


i=  1 s=  1 yj2n  I 2 


exp 


<Vis  ’ M - fcax.  - k/3y 


(41) 


From  (41)  it  is  easily  found  by  using  the  same  procedure  as  in  Section  4 that  the 
MLE  for  a and  /3  are  given  by  (16)  and  (17)  with  ^ replaced  by: 


M 


*i  " ™ *is  * 
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and  the  rms  errors  in  the  estimates  are  obtained  from  (19)  and  (20)  by  replacing 
2 . 
by 


2 

°t>  + 


n 

— 5 

A TB 


(42) 


2 

From  (42)  and  (19-20)  is  found  the  very  satisfactory  result  that  if  is  negligible, 
the  angular  accuracy  is  inversely  proportional  to  the  time  bandwidth  product  and  the 
ratio  between  signal  power  and  noise  energy.  This  result  can  be  put  in  an  even  more 
expressive  form:  If  Nq  is  the  noise  power  per  cycle  and  E is  the  energy  of  the  pulse 
of  length  T,  the  rms  error  in  the  angular  estimation  can  be  written  in  the  following 
form: 


a 


a 


2 


h2  * ^kh  i*  - v • 

' ' x y xy 


and  thus  the  contribution  to  the  rms  error  due  to  thermal  noise  (under  the  assump- 
tion of  high  signal  noise  ratio)  is  inversely  proportional  to  the  total  energy  in  the 
pulse  divided  by  the  noise  power  per  cycle. 


7.  CONCLUDING  REMARKS 

The  ultimate  accuracy  provided  by  a linear  or  planar  array  of  identical  sensors 
in  the  determination  of  the  direction  of  incidence  of  a plane  wave  has  been  studied 
for  the  ideal  case  of  the  absence  of  angular  interference.  If  the  phase  measure- 
ments are  affected  by  normally  distributed  errors,  the  optimum  processor  consists 
of  a linear  combination  of  the  phases  at  the  element  output  ports  with  weights  de- 
pending upon  the  array  geometry.  It  is  shown  that  conventional  multiple  interfer- 
ometer systems  are  substantially  less  accurate  than  the  "optimum"  processors 
introduced  in  this  paper.  The  conclusion  holds  also  in  the  presence  of  thermal 
noise. 

rtl enough  the  proposed  scheme  is  not  the  theoretical  optimum  when  terrain  re- 
flection is  present,  a very  limited  investigation  has  been  conducted  on  the  per- 
formance of  a small  array  of  sparse  elements  under  that  condition.  For  the  parti- 
cular case  studied,  the  maximum  errors  are  smaller  than  for  multiple  interferom- 
eter systems  and,  in  most  directions  of  the  limited  angular  range  investigated,  the 
errors  are  smaller.  How  the  scheme  compares  with  other  techniques  in  the 
presence  of  reflection  and  whether  in  this  situation  it  is  consistently  superior  to 
multiple  Interferometers  is  an  open  area  of  investigation  of  great  practical  interest. 
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Form  of  Estimator  for  Normal  Probability  Dansity  of 

Phase  Errors 


From  (15)  the  estimate  for  fi  is: 

N 

M = I (*,-k ax.-yp  yt) 

which,  introduced  into  (13-14)  yields,  by  recalling  (1)  and  (2): 

N N N N N 

T.  ^ x = Z X.  I Z (I  +koZ  (x.  - x)  X + k0  Z (y  - y)  X , 

i«=l  1 1 k-1  K N i=l  ‘ i=l  1 1 i=l  1 1 

N N N N N 

Z ^;X  =Z  y iz  16,  + kaZ  (x.  -x)y+k£Z  (y.  - >)  y,  . 

i=l  1 * k=l  K w i=l  1 i=l  1 ‘ i=l  1 1 


Notice  that  it  is  evidently: 

N N 

Z (x.  - x)  = Z (y.  - y)  = 0 . 
i=l  1 i=l  1 


On  the  basis  of  that  we  can  recogniie  that  (A2-A3)  takes  the  form: 


r M 58  T M 58 


$,ix.  - x)  = ka  M + k0  M , 
i i x ^ xy 

A.  A 

ip. (x.  - x)  = ka  M +k/3  M , 
ri  i xy  H y 


from  which  the  expressions  (16-17)  are  immediately  established. 


One  finds  the  parallel  result: 


Ef/2]  . 0 . 


Equations  (B3)  and  (B4)  indicate  that  the  estimates  of  the  angular  coordinates  a,  0 
are  unbiased. 


7-  **  ■ 


Appendix  B 


Mean  Value  and  Variance  of  the  Estimates 


Taking  the  average  value  of  (16)  and  recalling  (6-7)  one  gets: 

N 

L [/i  + ka  *i  + k 0yt]  - x)  Mj£  - (yt - y)  M ] 

E[a  ] = 1 i Zl 

k MM-  M* 


M M - M 
x y xy 


or,  through  some  manipulations,  recalling  (A4): 


. jt(m  m - m?  ) + fl  [m  m -mm 

1 ' y x xy'  ^ I y xy  y s 


MM-  M 
x y xy 


E[  a]  = or 
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The  variances  of  a and  0 provide  the  rms  error  in  the  angular  estimation.  By 
taking  into  account  (9)  we  have 

E [(<£  - Of)2] 

N N N 


E (x.  - x)2  M2  + E (y.  - y)2  M2  - 2 M M E (x  - x)  (y  - y) 
2 i=i  1 y i=i  1 ^y  x y i=i  1 


that  is  expression  (19).  Similarly  expression  (20)  is  obtained.  The  covariance  of 
the  estimations  of  a and  0,  entering  in  the  expression  of  the  joint  density  (22),  is 
found  to  be: 

E [(a  - a)  (0  - J)] 


N N N 


that  is  expression  (21). 


Appendix  C 

A Cramer-Rao  Bound  for  Unbiased  Estimates  of 
Direction  of  incidence 


In  this  appendix  we  address  ourselves  to  the  following  question.  What  is  the 
maximum  possible  accuracy  of  target  angular  location  that  can  be  obtained  by  pro- 
cessing the  observables  We  will  consider  only  unbiased  estimators,  that  is, 
those  not  having  systematic  errors.  We  restrict  ourselves  only  to  the  ideal  case 
of  no  angular  interference  (terrain  reflection,  and  so  forth)  and  no  thermal  noise. 

No  great  diffi  :ulty  would  be  encountered  however  in  taking  into  account  the  latter 
along  the  lines  of  Section  6.  No  hypothesis  on  the  form  of  the  phase  error  probabil- 
ity densities  (8)  is  made.  We  assume,  however,  the  phase  errors  to  be  uncorre- 
lated and  the  probability  density  function  to  have  everywhere  a second  derivative. 

To  simplify  notation  rename  the  unknown  parameters  as  follows: 


Cj  = a ; c2  = 0 ; c 3 = M . (Cl) 

With  this  notation  we  can  introduce  simply  the  so-called  Fisher  information  ma’.rix 
whose  elements  are  obtained  as  follows.  Take  the  opposites  of  the  averages  of 
the  second  derivatives  (with  respect  to  the  unknown  parameters)  of  the  logarithm  of 

g 

the  conditional  probabilities  of  the  observables: 


■1 

| 

1 

J 
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The  averages  in  (C2)  are  taken  over  the  observable  quantities  ar.d  the  c.  are  the 
"true"  values  of  the  parameters.  However,  it  will  be  apparent  in  the  sequel  that 
their  knowledge  is  not  necessary  to  evaluate  the  quantities  in  (C2).  According  to  the 
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Cramer- Rao  bound  therms  errors  ty.  in  the  unbiased  estimate  of  the  parameter  c. 
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satisfy  the  inequality: 

a?  > cofactor  of  J.j  . (C3) 


Our  task  is  thus  evaluating  the  set  of  J.^  given  by  (Cl). 

To  begin  with  notice  that  no  matter  what  the  explicit  form  of  p(6.)  in  (8),  we 
have,  from  (10): 
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Also,  since,  from  (6): 


(C4) 
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the  second  derivatives  of  do  not  contain  the  observables  ^ , Thus  the  second  term 


in  (C4)  is: 
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Thus  (C4)  is  simply: 
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The  Fisher  information  matrix  is  thus  explicitly  written 
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By  using  (C3),  through  some  simple  manipulations  we  obtain 

na2>  Je  f-  42l°-§P(e)  j 1 — ^ — g—  . 
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where  it  is  understood  that  the  probability  densities  are  conditioned  to  the  true 
values  of  the  unknown  parameter.  Equations  (C8j  and  (C9)  establish  the  maximum 
theoretically  possible  accuracy  for  an  angular  estimator,  when  the  phases  of  the 
sensor  measured  voltages  are  subject  to  instrumental  errors  6j.  whose  probability 
density  is  p(6.).  We  recall  that  this  result  is  valid  for  a single  target  (no  interfer- 
ence or  multipath),  assuming  identical  patterns  for  all  elements  (not  necessarily 
isotropic). 

For  normal  error  densities: 

i r n l \ - 1 


I E . d2  log  p(6)  ) 
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By  comparison  of  (19-20)  with  (C8-C9)  we  reach  the  conclusion  that  for  normal  den- 
sities the  phase  slope  reconstruction  algorithm  discussed  in  Section  3 is,  in  fact, 
the  best  possible  processing  scheme  in  the  sense  of  minimizing,  under  our  assump- 
tions, the  rms  error  of  the  estimate  of  the  angular  position  of  the  target. 
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Appendix  D 

Flimination  of  Phase  Ambiguity 

In  this  brief  discussion  on  the  elimination  of  phase  ambiguities  we  will  mainly 
refer  to  the  linear  array  discussed  in  Section  4. 

Let  us  denote  by 


A^  = (i  = 2,  3 ...  N) 


(Dl) 


the  observed  phase  difference  between  the  element  labled  i and  the  first  element. 
Consider  now  the  quantities: 


q.  = A 


x-  * x, 
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(D2) 


Because  of  the  very  definition  of  A it  is 


E[q.]  = 0 . 


(D3) 


where  the  averages  are  taken  with  respect  to  the  phase  measurement  errors.  Also 
it  is  not  difficult  to  establish  that  the  variances  of  the  quantities  (D3)  are: 
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For  i - 2,  (Dl)  gives: 


A *2  = *2  ' *1  ' 


(D5) 


and  we  assume  that  the  elements  2 and  1 are  close  enough  to  eliminate  the  possibil- 
ity of  ambiguity  in  the  phase  measurements,  although,  of  course,  measurement 
errors  will  be  present.  The  value  (D5)  will  belong  to  the  interval  -it,  it.  In  the 
array  considered  in  Section  4: 


for  all  i's.  Thus  we  have 


qi  = 


A ip.  - 2 A ip.  , 

Ti  Ti-1 


(i  = 3...  N) 


(D7) 


Once  (D5)  has  been  determined  in  an  unambiguous  way,  no  ambiguity  is  involved  in 
the  determination  of  to  be  used  in  (26),  if,  among  its  possible  values  differing  by 
multiples  of  2ir,  we  choose  that  for  which 


^3*2^2 
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Once  is  determined,  we  may  determine  by  repeating  the  procedure,  that  is, 
by  requiring  that  for  every  i: 

| A y.  - 2 A <,._1  | < n . (D9) 


This  leads  in  turn  to  the  specification  of  the  maximum  phase-measurement  error 
permissible.  To  show  that,  consider  that  from  (D4)  and  (D6)  we  get,  for  our 
geometry: 

Elqj2]  = 4 cr62  . (DIO) 

Assume,  that  with  probability  close  to  one.  the  quantity  q.  never  becomes  greater 
than  three  times  its  rms  value: 

A b.  • 2 A j s69{,  (Dll) 
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